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[Text] Abstract 


Research on an RF superconducting cavity at Beijing began in 1988. After three 
years of hard work, substantial progress has been obtained, and has been well 
recognized by experts worldwide. The preparation and procedures for the low- 
temperature superconducting experiment are described in detail. The paper also 
summarizes progress and accomplishments in RF (radio-frequency) superconduct- 
ing experimental technique, microwave and phase-locked measurement, post- 
purification of cavity, computer simulation design and computer control, and 
data acquisition and processing. Presently, an accelerating gradient of 8.6 
MV/m (Qo = 6.5 x 108) has been obtained in a 1.5 GHz niobium cavity operating 
in the CW mode. At 2 K, it was measured that Q, = 8 x 10%. These experimental 
results are at the cutting edge in the international field. 


I. Introduction 


RF superconducting technology has been used in accelerators for some time. In 
1965, a lead-plated superconducting cavity was developed at Stanford 
University to accelerate electrons. After more than two decades of effort, 
substantial progress has been made to improve the shape, material performance 
and accelerating gradient of superconducting cavities.! Presently, more than 
20 laboratories worldwide are using or building high-energy superconducting 
accelerators for heavy ions and electrons. 


The characteristics of a superconducting cavity include high gradient, large 
beam diameter, capability of operating in the CW mode, and low operating 














frequency.* There are advantages im using a superconducting cavity in a high- 
energy intense electron beam (e-beam) accelerator. A _ superconducting 
accelerator is particularly appropriate for a free electron laser (FEL) 
because it has very stringent requirements for beam quality, intensity, and 
mean power. In the late 80's, a scheme for using a superconducting cavity and 
an optical cathode high-intensity injector was presented with the objective of 
increasing the brightness of an FEL or TEV linear collider. ~ 


Today's superconducting cavity has excellent RF characteristics. In _ the 
laboratory, a single cavity has an accelerating gradient of 25 MV/m (Q, = 10%). 
Based on BCS theory, the accelerating gradient of a pure niobium (Nb) cavity 
can be as high as 50 MV/m.’ The experimental value is still far less than that 
value. Through a large number of experiments and continuous investigation, it 
is recognized that the surface characteristics of the superconducting cavity 
are critical to its overall performance. Processes such as multiple e-beam 
smelting, increasing RRR (residual resistance ratio), chemical washing, high- 
temperature treatment in ultrahigh vacuum, and artificial oxide-—film formation 
can effectively improve its performance. Nevertheless, field-effect emission 
is still the limiting factor preventing the acceleration gradient from rising 
further. Some large laboratories (e.g., KEK, DESY, CERN, and Cornell) are 
looking for new methods and techniques to modify the cavity surface in order 
to find a breakthrough in acceleration gradient.‘*:°°:’ 


The heavy-ion superconducting 
cavity group at Beijing Univer- 
sity has done a great deal of 
work on post-treatment of the 
cavity to raise the acceleration 
gradient. As a result, a series 
of effective measures has been 
established. Substantial prog- 
ress has been made in a rela- 
tively short period of time. We 
have measured a gradient of 8.6 
MV/m (in CW mode). At 2 K, a 
cavity quality factor of 8 x 10° 
has been obtained. This paper 
presents methods and techniques 
used to improve [the quality 
factor] Q and the acceleration Figure 1. Elliptical Cavity 
gradient Ege. 
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II. Design Analysis and PC-Based Program 
1. Design Analysis 


When an Nb cavity is running in a superconducting state, its resistance (Rg) 
is of the order of nQ. It is approximately 5 orders of magnitude lower than 
the R, of a copper cavity at liquid-nitrogen temperature.’ This extremely low 
resistance is a unique feature for Nb. An important point is that the shape of 
the superconducting cavity has very little impact on the RF power conversion 








efficiency. Hence, cavity design should address other issues such as effects 
due to dark are and higher-order modes (HOMs). Figure 1 shows an elliptical 
superconducting cavity. The beam tube has coupling openings and both ends have 
cutoff waveguides. A portion of the accelerating gradient has a transitional 
circular arc and it also has a large-aperture beam tube to effectively 
suppress the dark-arc effect and beam—induced field. The base mode, HOMs, and 
wake-field effects of a 1.5 GHz LE-type elliptical single harmonic cavity 
(f=1) have been analyzed using programs such as SUPERFISH, URMELT, and 1BCI. 
Relevant results have been published elsewhere.}°:!!'!* Figure 2 shows the wake 
field caused by a single focused pulse. An HP 8753C network analyzer was used 
to measure all the modes below 3.0 GHz associated with a 1.5 GHz Nb cavity at 
room temperature; the results are essentially consistent with the computer 
analysis. Table 1 is a comparison of measured and calculated values. 
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Figure 2. Wake Field Simulation of Elliptical Cavity 
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Table 1. Comparison Between Measurement and Computation 



































MODE URMELT Meas. MODE URMELT Meas. 

(MHz ) (MHz ) (MHz ) (MHz ) 
TMO10 1470.22 1467.70 T1110 2058.45 2043.54 
TMO11 2792.13 2800.40 TM111 2621.95 2585.17 
TMO020 2977.62 2975.84 TM112 2955.58 2974.58 
TE111 1889.80 1868. 32 




















2. Running Software Programs on PC 


In computer design analysis, an outstanding accomplishment was achieved. 
Programs such as URMELT and TBCI have been successfully transplanted onto a 
386-PC. To run a large program on a PC, the first issue is to solve the 
internal memory problem. The code and data segment of programs such as URMELT 
and TBCI are much larger than 640 K, which makes it impossible to run directly 
in a DOS environment. In the past, the conventional method has been to use 
overlays for the code and to open temporary files for the large data segment. 
However, this requires the careful analysis of each module and every 
executable statement in the source code. It takes an incredible amount of 
work. 


Our approach is to use virtual memory; this divides up a portion of the hard 
drive for virtual memory. When a 386-PC operates under a protected virtual 
address mode, it is capable of managing 16 M of real memory and 64 T (T=M x 
M) of virtual memory. This makes it possible to run large programs ona PC. In 
this case, the size of the program is determined by the capacity of the hard 
drive. On our HP-312 PC, with 10 M internal memory (2 M real and 8 M virtual), 
both URMELT and TBCI could be run normally. The speed is approximately one- 
third to one-half of that on the VAX. However, with multiuser time-sharing on 
a VAX, it is comparable in speed. Figure 2 shows some TBCI results obtained 
from the HP-312 PC. 


The fact that these programs can run on a 386-PC is extremely convenient for 
the research group, not only saving us time but also money. 


III. Cryostat and Ultrahigh-Vacuum System 
l. Cryostat 


The 1.5 GHz Nb superconducting cavity experiment has to be performed at 2~4.2 
K. To this end, a cryostat was constructed. After two years of effort, we 
completed the fabrication and installation of a system consisting of a large 
cryostat, a 100-liter liquid helium (He) Dewar, two 100-liter liquid-nitrogen 
(N.) Dewars, a liquid-He duct, an He recovery duct and gauges for temperature 
and liquid-He height. The construction of this apparatus ensures the smooth 
execution of superconductivity experiments. 


The 100-liter cryostat is the major equipment required to perform the 
experiment. The development of such a large low-temperature liquid-He device 
had never been reported domestically. Due to lack of information, the initial 
design was sent to DESY (Germany) for review by foreign experts. Our German 
colleagues came back with 10 major modifications. The design was then revised 
and made simpler and more practical. The new design stresses safety and the 
ultrahigh-vacuum requirement associated with the Nb superconducting cavity. 
Figure 3 is a sketch of the cryostat, which is approximately 175 cm high with 
a maximum outer diameter of 75 cm. The Nb cavity and the cryostat are on two 
separate vacuum systems. The cryostat uses a high-speed molecular pump and the 
system is required to provide 10°‘~10°° Pa of vacuum. The Nb superconducting 
cavity has an ultrahigh-vacuum system. 
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After about a year of operation, the 
entire cryostat system has essentially 
met the design requirements. Liquid He 
in the cryostat evaporates at a rate 
of 1 percent per hour. The 100-liter 
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liquid-—He Dewar evaporates at a rate 
of less than 1.5 percent per day. 
TO 10N nour = 

2. Ultrahigh-Vacuum System ; i 
The internal surface of the supercon- 
ducting cavity must be extremely SHI all 
clean. Otherwise, it has a significant ame Ss 5 
impact on Q and E,,,.. The 1.5 GHz Nb r4 | eee Sue 
cavity is vertically mounted in the ry GAIES- ens 
cryostat and is connected to the te tame TT a7 Slisue 
vacuum line in an overhead manner to r+ tap ao & 4 

: : ‘ 4 fa o- -) | 4 
prevent any debris from falling into entities 
the cavity. Both ends of the cavity 4 2 ea & 5 
are sealed with an indium (In) sealing ——— J | 74 | BEY 
surface. The inner diameter of the In ™ 44 e EB ox y Ly 
seal flange must be designed according 2 5 ~{{E= rad 
to the actual size of the beam pipe e gh Be 
flange. it must be highly accurate. 
The diameter of the sealing In wire is ~ 


4 
a mm. 


These two large-aperture sealed 
connections can directly influence the 
ultrahigh-vacuum achievable by the 
system. Especially when liquid He is 
in a superfluid state, these connections are the weak links. Two cryostatic 
absorption pumps and an ion pump provide an oil-free pumping system for the Nb 
cavity. The entire system is metal-sealed with ultrahigh-—vacuum metal valves. 
The RF feed—-throughs are 50-0 ceramic-seal coaxial connectors (total 99). An 
ultrahigh-vacuum jewel valve is installed at the end of the vacuum line as a 
control valve for helium treatment experiments. 














Figure 3. Liquid-He Cryostat 


This system can reach 10°°~10°’Pa of vacuum before using liquid N, to lower the 
temperature. After adding liquid N, and He, it can reach 10°°Pa. After repeated 
tests with liquid He, the system was found to be reliable for low-temperature 
physical experiments. The first phase of work was thus completed smoothly. 


3. Low-Pressure, Low-Temperature Experiment 


The Q value of a superconducting cavity varies with temperature.’ At 4.2 K, Q 
is (2~4) x 10° for the 1.5 GHz Nb cavity. At 2 K, Q = (6~9) x 10%. Hence, it 
is essential to achieve an ultralow temperature in order to reach a higher 
gradient. In order to reach an ultralow temperature, we employed a method to 
lower temperature and reduce pressure by a mechanical pump. After many tests, 
2 K was finally achieved in the 100-liter cryostat and Q, measured 8 x 10%. At 
2.12 K, liquid He is a superfluid, posing a serious challenge to the vacuum 
components in the liquid—He tank because superfluid He is extremely permeable. 
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After numerous low-temperature, low-pressure experiments, the vacuum inside 
the Nb cavity remained unchanged. This demonstrates the integrity of the In 
seals and other metal seals. It is rare in China for such a large liquid-He 
vessel to reach such a low temperature. 


4. Shielding of the Earth’s Magnetic Field 


The Earth’s magnetic field has a great deal of impact on the RF characteris-— 
tics of the superconducting cavity. Hence, it is necessary to establish a 
magnetic shield around the Nb cavity. We added several layers of thin film of 
a high-y material, which is an amorphous metal (Fe—-Ni-Si-B); the thin ribbon 
is 0.04 mm x 100 mm (thickness x width), total weight 6 kg. The axial and 
horizontal magnetic field inside the cryostat were measured by technicians 
from the Department of Geophysics at Beijing University using a model DM 2220 
Digital Magnetometer (U.S.-made). It was found that the axial magnetic field 
is one-seventh and the horizontal magnetic field is one-twentieth of the 
Earth’s magnetic field. This essentially satisfies the requirement for the 
superconducting experiment. 


IV. Post-Treatment of Nb Cavity 


After a superconducting Nb cavity is machined and welded, it needs to be 
polished mechanically, electrically, and chemically and then purified after 
putting on a titanium (Ti) film. As a result of these processes, a high- 
quality, high-thermal-—conductivity inner surface is produced. These processes 
help improve the quality factor Q and accelerating gradient E,... This involves 
materials science, chemistry and thermal treatment. In order to obtain the 
ideal RF superconducting conditions, the following work was done. 


1. Analysis of Nb Material 


The superconducting cavity inside an electron accelerator is usually made of 
high-purity thin Nb plate. Commercial-grade Nb cannot be used directly; 
several times of e-beam smelting are required to raise its purity. After 
raising purity, it is necessary to improve its thermal conductivity. The 
normal method is to evaporate a Ti coating on the Nb surface in order to 
remove the impurity in Nb by way of the combining process. 


We analyzed and compared the high-RRR Nb supplied by DESY to the e-beam- 
smelted Nb processed by the Institute of Nonferrous Metals. Table 2 shows the 
semiquantitative chemical composition of the two specimens. Sample #1 is from 
DESY and is not very impressive from the standpoint of impurity. Figure 4, an 
electron spectrum of the specimen, shows that the surface is covered by a 
layer of Nb,0, several hundred angstroms in thickness. After repeated RF 
superconducting experiments, it was verified that this oxide layer could 
stabilize the RF characteristics of the cavity. Presently, this high-quality 
oxide film is being added to the purification process at KEK and CEBAF as 
well. 








Table 2. Semi-Quantitative Chemical Analysis of Two Nb Specimens 
Sample #1 [Figures in percent] 




































































Ni Zr Rh Sn Tl Ta Mg Pt Ti 
0.001 | 0.003~ | 0.01 0.001 0.001 0.1 0.001 0.002 0.001~ 
0.01 0.01 
Sample #2 [Figures in percent] 
| ue 2x n | a Pt cu rif 
| 0.01 0.003 0.001 0.003 0.001 0.001 0.003 | 
xoos«z e8-3-89 «(re Bt wot — 
2. Chemical Polishing REGION ) aan hd 
xPS AL CAE,SO S.m'57 396 SEC 0.050 STEP 
The 1.5 GHz Nb cavity manu- wt 





160] 


factured by Dornier Corpo- 
ration is made of high-pu- 
rity Nb with RRR = 300. 
Usually, the thermal break- 


down threshold is higher Cnt oe | 
when RRR is high. It will ~~ ,,,3 

be easier to attain a high-— RATE | 
er accelerating gradient <1I0? 4, 


Eacc; however, it is more 
difficult to chemically 60 unas 
clean this high-RRR Nb. In 

the acid washing process, rT { - 

it is easy to incorporate = 
H* or phosphate into the Nb 
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face characteristics. In 
order to avoid this situa- 
tion, we conducted numerous chemical polishing experiments and finally 
established an effective polishing method. 


Chemical polishing was done according to the method used at DESY and KEK.?° 
Figure 5 [photos not reproduced] shows scanning electron micrographs of a 
specimen before and after chemical polishing. It is obvious from the pictures 
that the crystals are well aligned on the Nb surface after chemical polishing. 
The polishing solution is composed of a fixed ratio of HF (40%), HNO, (65%), 
and H,P0O, (85%). The number of washes and the time of exposure must be 
carefully controlled. The Nb cavity needs to be cooled during the process. 
Cleaning after acid polishing is very critical. The resistivity of the 
deionized water used must be monitored frequently; especially in the final 
rinse, it should be at approximately 18 MQ. The cavity is mounted on a special 
support during acid polishing. It can be rotated to allow the acid to come 








into contact with the cavity wall evenly and to make the discharge of the 
solution easier. Acid-resistant flanges and valves are used on both sides of 
the beam pipe. The first chemical polishing prior to thermal treatment must 
remove a surface layer 20~30 uw in thickness. After thermal treatment, it needs 
to be chemically polished one more tine to remove another 5~10 yu 


3. High-Temperature Treatment 


We began conducting RF superconducting experiments in October 1990. In the 
first two experiments, only chemical polishing was done and the results were 
not satisfactory. After reviewing the data obtained in KEK, DESY, and Cornell, 
our preliminary conclusion was that stress in the cavity had to be relieved. 
stress comes from machining and welding. According to a 1989 Cornell 
publication, it was reported that the Q, for their 1.5 GHz Nb cavity was two 
orders of magnitude lower than the expected value and could not be improved 
until the stress was relieved by heat treatment. DESY also presented the same 
issue in a 1989 meeting. In reference to their experience, a thermal treatment 
procedure was formulated as follows: The heat-resistant is done in a vacuum 
furnace with a vacuum better than 1 x 10°°Pa and at a temperature of 850~900°C 
for one hour. The vacuum pumps should not be turned off during cooling. The Nb 
cavity is placed in a Ti box during the heat treatment. Impurities in Nb are 
evaporated from the surface at high temperature and may be readsorbed during 
ling Ti box, impurity ions are easily adsorbed by Ti to keep the Nb 
cavity from secondary contamination. 


cooling. With a 


The Nb cavity used in our third liquid-He-temperature experiment was thermally 
treated; this is the very cavity with which we had a major breakthrough. The 
measured Q, was of the order of 10° (at 4.2 K). With reduced pressure and lower 
temperature, Q, was found to be of the order of 10%. These numbers are 
consistent with results reported in other countries and they are also in 
agreement with theoretical values. Furthermore, we also noticed that the heat- 
treated cavity was very reproducible. We did not see Q, drop due to temperature 
excursion. ?° 


4. Ultraclean Installation 


In order to prevent contamination, the cleaned cavity must be installed in an 
ultraclean environment. After chemical polishing, the cavity is filled with 
ultrapure water to prevent any impurity from entering the cavity. It is then 
immediately sent to an ultraclean room for installation. Prior to installa- 
tion, the cavity is rinsed again with ultrapure water. The ultraclean room is 
equipped with an air curtain and air filters. It also has a Class-—100 work- 
bench for two to ensure that there are fewer than 100 dust particles more than 
0.5 wm in size per cubic foot over the workbench. The number of such dust 
particles in a cubic foot of normal air is 210,000. When installing the main 


coupler, pure nitrogen evaporated from the liquid-nitrogen Dewar is fed into 
the cavity through the heating tube to create a positive pressure inside the 


cavity. This is to prevent air from coming in contact with the cavity surface 


to minimize the chance of any impurities entering the cavity. After the cavity 
is connected to the main coupler, nitrogen is injected from the main coupler. 








These two components are then connected to the vacuum system. During this 
process, a positive pressure is always maintained inside the cavity. 


After turning on the ion pump, the entire system is baked for approximately 20 
hours, and will reach the desired vacuum in a week or so. 


V. Microwave and Phase-Locked Measurement Technique 


Microwave measurement is a key step in RF superconducting experiments. Under 
superconducting conditions, the resistance R, of the 1.5 GHz Nb cavity is less 
than a few hundred ohms, which is five orders of magnitude lower than the R, 
of a copper cavity. Hence, the quality factor of the superconducting cavity is 
as high as 10°~10°. The half-power bandwidth of its resonant curve is only a 
few Hertz wide and the field-decay time constant is 30~300 ms. These proper- 
ties make the measurement of RF parameters and search for the resonant point 
very difficult. From the very beginning, we recognized that microwave measure- 
ments involving such a high-Q superconducting cavity were never attempted in 
China and it would require the appropriate test methods and instrumentation 
with high accuracy and stability 


1. Instrument Setup 
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Figure 6. Microwave Phase-Locked Loop 





Figure 6 shows the microwave measurement circuit used in the low-temperature 
superconducting experiment, and essentially shows our instrument setup as 
well. The synthesized signal generator is an HP 8663A, which has a frequency 
stability of 10° ’°/day. The output signal has a resolution of 0.4 Hz at near 
1.5 GHz, which is less than the half-power bandwidth of the superconducting 
cavity. The signal generator is capable of scanning and phase-locked tuning. 








When searching for the superconducting point by frequency scanning, a signal 
generator must have the kind of high resolution and stability as those 
described above. 


The bandwidth of the four-channel digital oscilloscope, an HP 54503A, is DC- 
500 MHz. The horizontal time axis has 501 points. (The HP-IB remote data 
acquisition system has 1,024 points). The vertical voltage axis is eight-bit. 
Its accuracy is 1/2° = 0.4 percent. (The HP-IB remote data acquisition system 
is 16-bit with an accuracy of 1.5 x 10°°.) It can automatically detect the 
maximum and minimum signal. Furthermore, it is capable of measuring the 
voltage at a specific time by remote control. It also can determine the time 
at which the signal is at a specific voltage. The oscilloscope also has memory 
for storage 


The power meter is an HP 43/7B and its measurement range is: 100 kHz--50 GHz in 
frequency, -—/70 dBm to +44 dBm in power, and 0.1 percent in accuracy. The 
frequency counter is an HP 5350B. It has an accuracy of 1 Hz near 1.5 GHz and 
the gate is 200 ms wide. 


The HP PC-312 controller has two CPUs. It not only can function as a 386 
microcomputer but also can be connected to thie signal generator, oscilloscope, 
power meter, frequency counter, and laser printer through an JEEE-488.2 
interface to form a fully automated data acquisition, processing, storage, and 
print-out system 


The 20-W RF solid-state amplifier is the final-stage power amplifier of the 
measurement system. It was constructed by Institute 54 of the Ministry of 
Machine-Building and Electronics Industry (MMEI). The frequency range of the 
amplifier is 1.4 GHz~1.6 GHz. When Q, of the Nb cavity is of the order of 10%, 
the maximum power output meets the experimental requirement. 


The system employs directional couplers such as HP 778D (20 dB), Narda 4242-20 
and TT3-W (10 dB), and an HP 8470B crystal detector and a DMM-6-1500 dual 
balance mixer (DBM). 


2. Phase-Locked Loop (PLL) 


The half-—power bandwidth at the resonant point of the superconducting cavity 
is only a few Hz. However, the frequency drift caused by external mechanical 
vibration is approximately 20~30 Hz. The frequency drift during cooling is 
approximately 0.1 MHz/K. To this end, the microwave measurement system 
requires a PLL technique to allow the frequency of the signal generator HP 
8663A to track the frequency change of the cavity in order to ensure that the 
superconducting cavity remains in a resonant state. 


Figure 6 shows the PLL measurement circuit. The microwave signal generated by 
the HP 8663A (10 mW) travels through the 20-W power amplifier, isolator, and 
variable attenuator to reach the HP 778D directional coupler. At that point, 
the main signal is directly fed into the main coupler of the superconducting 
cavity from the exit of the directional coupler. In addition, a -20 dB signal 
is extracted from the incident wave to be transmitted to one port (LO) of the 
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DBM after passing through a mechanical phase shifter. Another microwave signal 
is extracted from the superconducting cavity to be transmitted to the RF end 
of the DBM. A phase-error signal will be sent after comparing the two signals. 
This error signal goes through a PI amplifier and is then transmitted to the 
FM port of the HP 8663A to form a stable feedback-control circuit. 


3. Main Coupler and Pick-Up Probe 


The microwave power is fed into the superconducting cavity through the main 
coupler. The pick-up probe extracts the microwave signal from the cavity for 
measurement. As far as power input and pick-up signal are concerned, the 
external quality factors, Q,,, are defined as follows: 


Qoxe, = 2NU/W, (input portion), Qoxe, = 2NU/W, (pick-up portion) (1) 


where U is the energy stored in the cavity, and Wg and W, represent the energy 
loss per cycle at the internal resistance of the signal generator and that in 
the signal pick-up circuit, respectively. It is known’ that the quality factor 
with load, Q,, is related to [the unloaded quality factor] Q, as follows: 


1/0, = 1/Q, + 1/Qexe, ’ 1/Qoxe, (2) 


In order to express the degree of coupling, i 
let us introduce the coupling factors for the COaxiaL CONNECTOR 
input circuit and pick-up circuit as f, and 


Bp. 





B, = 0o/Qoxt, + B, ° 0/Qoxt, (3) 


Substituting §£, and B, into equation (2), we 


have: MAIN COUPLER 








Q = Q, (1 + By + Bo) (4) 


The quality factor of the cavity with load, 
Q,, can be measured experimentally. The 
unloaded quality factor, Q, can only be PICK UP 
determined after 6, and f, are obtained. 
Figure 7 is a schematic diagram of the main 
coupler. The Q) value is approximately 10°~ 
10°. In order to feed as much microwave power 7 L' 
as possible into the cavity, it is required COAXIAL CABLE ~—— | 
that Q,,,. of the main coupler be comparable 
to Q . Since Q,,,. of the pick-up probe is two Figure /. Main Coupler and 
orders of magnitude higher than Q), the Pick-up 

magnitude of the signal picked up by the 

probe has little effect on the measurement of Q,). 6, of the probe can be 
neglected in calculating Q). 
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The external Q values are designed to be: 0, = (1~3) x 10%, Qoxe, = 1 x 10% 
and the length of the main coupler antenna and that of the probe are 
determined accordingly. These values were first estimated using a computer and 
then adjusted experimentally; finally, they were verified during superconduc-— 
tivity. If necessary, the devices will be reinstalled; this is a delicate 
piece of experimental work. At liquid-He temperature, the value of Q varies 
due to thermal expansion and contraction. The pick-up probe is made of a 50-0 
hard coaxial cable. During cooling, the insulating layer between the inner and 
outer conductor contracts several dozen times more than the metal. The result 
is that its Q,,,. is two orders of magnitude higher than the designed value at 


liquid—He temperature. We encountered similar problems with the main coupler. 


The antenna of the main coupler is a piece of 50-Q] coaxial cable. There is a 
vacuum between the inner and outer conductor. One end is connected to the 
power feed line through a 50-Q coaxial ceramic seal. This ceramic connector 
can only operate down to 20 K. The inner conductor is made of molybdenum. The 
antenna and the inner conductor are threaded together. Initially, the antenna 
was made of copper. After the experiment, it was found that the antenna almost 
came loose and there was only one thread left connected. This is because 
copper has a higher thermal expansion coefficient than molybdenum. Through 
experimentation, the design was modified and materials were changed, in order 
to make experimental values obtained at room temperature consistent with those 
obtained at low temperature. 














TO FOVER METER 
The electric field distribution, E;, sn asan C 
in the cavity was calculated using Hk 7780 } 
URMELT. Based on the relationship |5EFF °sct4ATor 











between probe voltage, electric field 
in the cavity and probe length, it is 
possible to estimate the length of the 
probe. Figure 8 shows the experimental 
method to measure the lengths of the 
antenna and the probe. An auxiliary 
antenna is set at one end of the cavi- 
ty and its length is adjusted based on 
the reflective scan signal to make 

8B = 1. The output power of the signal 
generator is measured with a power 

meter. Then, the pick-up power of the 7 Wy SENSO 
main antenna at the other end of the er 
cavity is measured. The length of the 
antenna is adjusted to obtain a suit- 
able power pick-up to output ratio. 
Thus, the desired external Q values can be obtained. In our experiments, 
Qexe, = 1x 10°, Ogre, = 2 x 10'°. The latter is one order of magnitude lower than 
the designed value because the phase comparator requires the probe to provide 
a signal higher than 5 dBm. 
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Figure 8. Principle of Q,, Measurement 
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VI. RF Superconducting Experiments and Results 


RF superconducting experiments began in October 1990. The first experiment was 
to test the entire low-temperature system. The low-temperature and vacuum 
system was found to perform well; it was safe and reliable to operate. The 
system is capable of pumping the Nb cavity and the cryostat jacket to 10°°Pa 
and 10°‘Pa, respectively, prior to cooling. After two days at liquid-N, 
temperature, the temperature of the Nb chamber reaches 160 K. During a 
superconducting experiment, liquid He is transported from the Dewar to the 
cryostat in a pipe until the entire Nb cavity is submerged. At the same time, 
the recovery pipe is turned on to receive He in the low-temperature workshop. 
Afterward, the liquid-He pipe is removed. The temperature of the Nb cavity is 
4.2K. The level of liquid He is above the upper flange of the beam pipe. The 
temperature and liquid-He level are monitored with an Fe-Rh thermistor and a 
liquid—He-level meter, respectively. A 15-liter mechanical pump is attached to 
the cryostat with a three-way connector. Its exhaust is connected to the 
return. The reduced-pressure, low-temperature experiment can proceed by 
turning on the mechanical pump and controlling the valve. The lowest tempera- 
ture achievable is 2 K. 


1. RF Superconducting Results 
10° 
The RF superconducting experiment is designed to \ 
investigate the RF characteristics of the 1.5 
GHz Nb cavity in the superconducting state. This 
kind of L-band resonant cavity is used in a 
linear accelerator. Experimentally, the unloaded 
Q value was measured as a function of accelera- 
tion gradient, E,,,, and temperature, T. Tempera- 
ture was varied by cooling under reduced pres- 
sure. Q values of the cavity at 6.5 K, at 4.2 K, 
and at several points between 4.2 K and 2 K were 
measured. Such measurements were made over and 
over again with excellent’ reproducibility. [ 
{ 
2 


ie 





Figure 9 shows a curve plotting surface resis-— 

tance R, against T,/T. The relationship between 10" 
R, and Q, is G = Rg x Q, where G is the geomet- 

ric constant of the cavity. For this cavity, T/T 

G = 260 Q. T, is the critical temperature and is Figure 9. Curve of R, Vs. 
9.2 K for Nb. The solid line in Figure 9 is the T,/T 
expected R, versus temperature curve derived from 

BCS theory. The triangles are experimental points. The theoretical values are 
in very good agreement with experimental points in most areas. At very low 
temperatures, the experimental points deviate from the theoretic curve and 
gradually approach a fixed value. This is normal and the value is known as the 
residual resistance, Rrpe., of the superconducting cavity. This is also the 
difference between RF superconductivity and dc superconductivity. The 
magnitude of Rre, is dependent on Nb surface characteristics.’ 





Figure 10 shows the Q, versus Eg, curve. A power experiment was performed at 
4.2 K and 2.25 K. Q and E,., were measured at various input power. As shown in 
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the figure, Q, decreases with in- Q, 


creasing Eg.,~. The rate of decline io 
and the presence of an obvious ae 
turning point are criteria to judge 


the RF characteristics of a cavity. 
The Q) versus Eg, curve should be 
relatively flat for a good cavity. 
At a higher gradient of acceler- 
ation, Q, should drop or show an 
apparent turning point. No turning 
point was found in our experiment. 0 ' 2 ; 4 $ ; . , . 
However, Q, did fall somewhat. In Ee (MY/m) 

our initial power experiment, a Figure 10. Q Vs. Eacc 
gradient of acceleration of 3 MV/m 

was obtained. After repeatedly 

improving techniques such as chemical polishing, thermal treatment, reduced 
pressure cooling, and microwave PLL, it has been increased to 8.6 MV/m. This 
value is already at a world-class level. Currently, we are preparing to 
conduct the more difficult He experiment to further improve the Q) versus Eg, 
curve to obtain a larger Q) at higher Eg,,<. 





) 
10° } 








In the figure, E,.. is calculated from the following: 





Eyoc = (¥2/0.1)/102 x Q, x (4B/(1+B)*) x P, (5) 


where 0.1 is the effective acceleration length of the cavity (m), (102 x Q) 
is the impedance of the cavity (M2), P; is the incident RF power (W) and E,,, is 
in units of (V/m). Values of # and Q) are discussed in detail in the following 
section. 


2. Microwave Measurement Technique Under Superconducting Conditions 


As mentioned earlier, the unloaded quality factor, Q), of the 1.5 GHz Nb cavity 
is 10°~10° and the half-—power bandwidth is only a few Hz. The field—decay time 
constant, however, is often as long as 30~300 ms. This makes the measurement 
of RF parameters and location of the resonant point very difficult. Further- 
more, because our experimental cavity does not have a frequency tuning 
mechanism, it makes it more difficult to determine the resonant frequency. 


(1) Locating the Resonant Point by Frequency Scanning 


Because the HP 8663A synthesized signal generator has a very high frequency 
stability and a 0.4 Hz frequency scan resolution, precision scanning is 
carried out in every 100 Hz region between 1.46 GHz and 1.47 GHz at a rate of 
100 ms/Hz. After repeated precision frequency scanning, the superconducting 
resonant frequency was captured. After this frequency is obtained, the 
generator is switched from scanning to CW mode. In addition, the PLL begins to 
operate. Under the control of the PLL, the superconducting cavity operates 
steadily in a resonant state. 
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Figure 11. 8 and Reflected Pulse 
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(2) Measurement of 8 by Reflection 


The value of 8 is measured by using the transient response of the resonance 
circuit. When a pulse-modulated RF signal is fed into the cavity, the electric 
field inside the cavity undergoes a build-up and a decay process. By 
subtracting the transient incident voltage pulse from the transient cavity 
voltage signal, we have the transient reflected voltage signal. The shape of 
the reflected voltage pulse varies with 8. Figure 11 shows three different 
reflected pulses corresponding to 8 <1, 8 = 1, and 8 > 1. From transmission 
line theory,*® we have 


B = X/(2 — X) (6) 


where X = b/a; b is the height of the latter peak in the reflected pulse and 
a is the height of the earlier peak of the reflected envelope. 


Ae ewelting trigger IHIP of Peking University 
The transient re- . 26 Mer 1991. 15:28:56 
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close to 1 as pos- 

sible to minimize 

any measurement error. The best way is to use a main coupler with an adjust- 
able antenna. Figure 12 shows the reflected pulse during a superconducting 
cavity experiment. 
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(3) Measurement of Q by Decay 


Experimentally, the loaded quality factor Q, can be measured. With the 
corresponding §, Q, can be calculated. Q,; may be measured by frequency scanning 
or decay. When Q is as high as 10°~10°, f, can no longer be accurately measured 
by scanning. However, the decay technique is suitable for measuring a high Q 
value. When a stable electromagnetic field is established in the cavity, the 
RF signal is cut off to measure the decay curve of the field. A time constant 
r can be determined from the decay curve. The value of Q, can be calculated 
from the following: 


Q, =x" xf)xT (7) 
where f, is the resonant frequency. 


The other curve in Figure 12 is the free decay signal of the electromagnetic 
field in the cavity. 


3. Control of Experimental Apparatus and Data Acquisition 


In terms of microwave measurement technology, we have successfully located the 
resonant frequency of the high-Q cavity and kept the superconducting cavity 
operating in resonance with stability using a PLL technique. These methods are 
new to the microwave measurement field in China. In the meantime, micro- 
computer control and automatic data processing were made possible for the 
system shown in Figure 6. Programs were written for an HP PC-312 microcomputer 
to control the entire measurement system (signal source, oscilloscope, 
frequency counter, and power meter) to make measurement accurate, intuitive, 
fast, and convenient. 


The HP PC-312 has an HP BASIC controller which has its own CPU on board. The 
controller board has a 2.5 MB BASIC dedicated memory for remote control and 
data acquisition through an HP-IB interface. The program written is specifi- 
cally for the Q values of the RF system and its capabilities include on-line 
calculation, screen (result) printing, data storage, data display, and post 
data processing. A menu is used to select and modify parameters. It is simple 
to operate. The software is comprised of top-down modules. It has a solid 
structure and has a great deal of fault-tolerance. Its user interface is also 
friendly. The program can perform the following: 


(1) Measure Q by way of frequency scanning or decay and display results on the 
oscilloscope screen. 


(2) Print any position of the oscilloscope screen with a laser printer in four 
different sizes. Several figures may be printed on the same page. 


(3) Oscilloscope waveform data is stored on disk for future recall, recalcula- 
tion, and post treatment. 


(4) Redisplay and computation of stored waveform. 
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Figure 13. The Function of Program Control 


CURRENT STATUS 





Figure 13 shows a 
brief outline of the Mode Reflect Pick-u niin 
function of the con- rHaNNEL >| mM 

trol program. Figure 

14 shows the main 


wee of the progran. Main MENU for Measurement of 0 
Figure 15 is a pic- wt--Chenee 
vale an 
ture of the 1.5 GHz CHANNEL-- Change 
Nb superconducting Soka’ heim 
cavity. Figure 16 PRINTER--Output 
SAVE--DATA to File 
shows the assembled SHOW----- Old-DATA 
Nb cavity. Figures 


17 and 18 are photo- 
graphs of the micro- 
wave PLL measurement 
system and the cryo- 
stat, respectively. 


The upper part of User 1. Caps Command 
calculated Q and 8. 
Qo, can be derived 
based on these num-— 
bers. Program controlled automatic measurement is very convenient and timely. 
It improved our efficiency and success rate. 





Figure 14. Main Menu of Program 


VII. Conclusions 


The superconducting cavity group at the Institute of Heavy Ion Physics of 
Beijing University has made significant progress in the study of the RF 
superconducting cavity in a relatively short period of time. This research 
will serve as a basis and technical foundation for future use of the RF 
superconducting cavity and development of superconducting high-brightness 
electron injectors. 
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[Text] Abstract 


The electron beam (e-beam) from the wiggler, after passing through the 180° 
bend, enters the energy and time spectrum diagnostic system. The 180° bend is 
designed for achromatic and isochronous transport of the e-beam. Isochronous 
transport satisfies the diagnosis of the time resolution of micropulses and 
achromatic transport satisfies the diagnosis of the time resolution of macro- 
pulses. The time-dependent macropulse energy spectrometer is comprised of four 
matching lenses, a two-slit device, a sweep deflector coil, a 90° symmetric 
bifocal analyzing magnet, a fluorescent target at the focal plane, a video 
camera system with an image intensifier, and an on-line computer-based 
graphics processing system. Its energy resolution is better than 1,000 and 
time resolution is higher than 60 ns. 


[Introduction] 


Figure 1 shows the layout of the diagnostics for the time-dependent energy 
spectra of micro and macropulses of the e-beam of the Beijing FEL. Limited by 
available space, the e-beam from the wiggler must be bent by 180° before 
entering the time-dependent energy spectrometer system. This paper describes 
the 180° bend and the macropulse diagnostic system 
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Figure 1. Beamline Configuration for E-Beam Diagnostics of Beijing FEL 


I. Design of the Isochronous and Achromatic Bend 


The 180° bend may be designed for either isochronous or achromatic transport, 
or both. Isochronous transport satisfies the diagnostics of the micropulse 
time spectrum I(t) and achromatic transport satisfies the diagnostics of the 
time-dependent macropulse energy spectrum I(E,t). 


For a system comprised of a horizontal bending magnet, a quadrupole lens and 


a drift segment, its first-order six-dimensional transport matrix is as 


follows’: 


M, M;, 0 0 0 My! 
M,, M22 0 0 0 M2. 
0 0 M, M, 0 0 (1) 
0 0 M, M, 0 0 
Meo, Msp 9 O 1 Mog 
fo 0 0 0 0 1) 
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Furthermore, when it satisfies the following conditions: 


oO 


M,, = 
te s- (2) 
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it is achromatic. This means that the emissivity of the e-beam does not 
increase as it passes through the system despite the fact that its energy is 
dispersed. When 


Ms, = 


Ms. = 
Ms, = 


(3) 


oo oO 


the system is an isochronous transport system. This means that the deviations 
in trajectory and energy of the e-beam are independent of the initial motion 
on the bending plane and the longitudinal length of the electron cluster 
remains unchanged. 


In engineering terms, equations (2) and (3) imply that the system must have at 
least five optical elements, such as quadrupole lenses, to provide complicated 
adjustments and diagnostics to put the system in isochronous and achromatic 
operation. Obviously, this is very inconvenient. To simplify the desig. and 
for ease of use, the 180° bend was designed symmetrically. Equations (2) and 
(3) can be simplified and the transport matrix elements at the center of the 
system (on the symmetry plane) can be made to satisfy the following conditions 
(proof for isochronous and achromatic transport is omitted). 


M,, = 0 (4) 


M,.M., - M,,M., = 0 (5) 


Mz, + Mz¢(Ms,M2, - Mc2M,,) = 0 
In practice, two optical elements for adjusting the beam parameters can be 
eliminated. Considering the fact that the streak camera used to measure the 
Cerenkov radiation of the micropulse is a block-type measurement, i.e., 
micropulse and macropulse diagnostics cannot be done simultaneously, the 180° 
bend was designed to have two pairs of symmetric quadrapole lenses (each pair 
in series). When performing macropulse diagnostics, only one pair is required. 


This kind of symmetric bend has several advantages. For example, when the 
system operates in the achromatic mode, equation (4) becomes: 


(Hy2Q2, + Hay(Q,, + 1202,)] (Aye + Hog dy) (6) 
+ [(Q,, + 1,022) Ha, + Q22H22] Hog + Hz, = 0 


where H, W, and Q are the transport matrixes of the first bending magnet, 
second bending magnet, and quadrapole lens, respectively; 1, and l, are the 
distance from the exit of the first bending magnet to the inlet of the lens 
and the distance from the exit of the lens to the inlet of the second bending 
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magnet, respectively. By substituting the parameters of or system into 
equation (6), we have 


24.787 ksin kl — 0.03965 Suk —2.6895cos kl - 0.707 = 0 


where k = (G/p)}/2, G is the gradient of the lens (Gs/cm), p is the electron 
magnetic rigidity (Gs-cm), and 1 is the effective length of the lens (1 = 10 
cm). A unique solution of k may be obtained from this transcendental equa- 
tion. In practice, the excitation current of the lens, Ig, is set based on the 
magnitude of the excitation current of the bending magnet (i.e., electron 
magnetic rigidity), Ig, to make k = ky which puts the system in the achromatic 
mode. Furthermore, the Ig~Ig curve also serves as an operational curve which 
is precalculated. When the system is set for the isochronous mode, there is 
another set of Ig,, Ig2~Ig curves. Similarly, these curves can be calculated in 
advance. 


In mid-1990, the Beijing FEL group requested that the micropulse measurement 
system be installed behind the first 90° bend. The original achromatic design 
for macropulse diagnostics was not affected by this decision. However, the 
symmetric isochronous transport design would not meet the requirements for 
micropulse diagnostics. In operation, careful adjustments of the parameters of 
the two lenses in the first 90° bend are required to achieve isochronous 
transport in the 90° bend. Computation done with the TRANSPORT program! 
indicates that isochronous transport is still possible in this case to meet 
the requirement for micropulse diagnostics. Furthermore, the calculated 
results can provide guidance for operation. Figure 2 shows the calculated 
e-beam envelope in achromatic and isochronous transport. The initial state of 
the e-beam is that it is in the center of the wiggler. The e-beam has two 
waists and the diameter of the waist is 1 mm. The normalized emissivity is 
30x mm—mrad. The focusing in the y direction provided by the B, component of 
the wiggler field is linearized as follows: 


y cos (z/ .) B .sin(z/B.)|[ y | 
y’ | (1/9 .)einte/B cos (2/6 .) y’ |oco 
where Be = N dy/2n. 


N is the number of periods of the wiggler field and Ay = 3 cm. 
II. Macropulse Time-Dependent Energy Spectrometer 


The system consists of four matching quadrupoles (after a 180° bend), a double 
slit device, a deflector sweep coil, a 90° symmetric bifocal analyzing magnet, 
a fluorescent target at the focal plane, a video camera with image intens- 
ifier, and an on-line computer graphics system (Figure 1). The following is a 
description of the design of major components. 
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1. 90° Bifocal Analyzing Magnet 


After achromatic bending of 180°, the e-beam is transformed into a double 
waist (object) at the horizontal slit of the double-slit device. An image is 
formed on the fluorescent target after passing through the 90° bifocal 
analyzing magnet. This analyzing magnet is designed to be symmetrically 
bifocal for more stable operation. Its bending radius R = 400 mm, horizontal 
dispersion Dy = 4, and amplification factor M = 1. With a full object slit 
width of 1.6 mm, the energy resolution is 1,000. This means that 10 measure- 
ment points can be obtained on the focal plane for an e-beam with a +0.5% 
energy dispersion. In this case, the full image width is 16 mm. The designed 
maximum e-beam energy dispersion is +3%, which corresponds to a full image 
width of 96 mn. 
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Figure 2. Envelope of Electron Beam 


In order to ensure magnetic field uniformity in the working area (better than 
10°°), the width of the pole is four times the air gap between poles. The 
operating magnetic flux is 0.1~0.35 T. In order to ensure the same effective 
hard magnetic boundary at different field strengths, the boundary is designed 
to be a Rogowski approximation boundary with magnetic shield clamping. The 
incident and emergent angles of the analyzing magnet are designed to be 
tunable over the range of 25° + 15° so that the actual incident and emergent 
angles may be adjusted according to the measured edge field distribution. 
Furthermore, this arrangement will make it more convenient in the future to 
compensate for aberrations near the poles with a circular arc boundary 
provided by a hexapole field. Figure 3 shows the uniform field distribution 
calculated using the POISSON program.* The data (Table 1) indicates that the 
uniform region is over 4 cm wide. 
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Computation done with the 
TRANSPORT program shows 





that when the normalized 
emissivity of the e-beam is 
30x mm-mrad, the object 
slit can vary between 1.2 
and 1.6 mm without affect-— 
ing the transport charac-— 



































teristics of the system. 
The corresponding resolu- 
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meter is 1,300~1,000. 


Figure 3. Calculated Results of the 


POISSON Program 























Table 1. 

ij x (mm) —30 -15 +15 +30 

| AB/Bo 5.4 x 10% 1.6 x 10“ 5.7 x 10° 5.8 x 10 
y 

2. Sweep Deflector Coil y 


The diagnostics of the time-dependent 
macropulse energy spectrum I(E,t) are 
accomplished by a vertically sweeping e- 
beam generated by a sweep deflector coil 
within a macropulse. In this case, the 
vertical beam size on the planar fluo- 
rescent target represents the width of 
the macropulse. The macropulse width of 
the Beijing FEL is 2~5 ws. Hence, the 
sweeping magnetic field rises very rap- 
idly. In order to avoid a very large 
inverse peak voltage between the magnet- 
ic pulsed power supply and the excita- 
tion coil, a low-inductance coreless 
saddle-shaped deflector coil is used to 
generate the uniform bipolar magnetic 
field. The coil, whose structure is 
shown in Figure 4, is split into two 
symmetric halves which are_ tightly 
wrapped around the beam pipe. The angu- 
lar distribution of the number of wind- 
ings varies in a manner close to a co- 
sine function. Since the coil inductance 
increases with the square of the number 
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Figure 4. Scheme of Slow-Deflector 
Sweep Coil 


of windings, the number of windings should be as small as possible. However, 
the fewer the number of windings, the poorer the uniformity of the bipolar 
field and the larger the error of the diagnostics of the time-dependent 
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spectrum. We were required to use the smallest number of windings provided 
that the uniformity of the field is less than 1 percent. A variety of schemes 
were compared by computation using the MPVP program’ and by actual measurements 
using a simulator. The final design chosen has a total of 52 windings and its 
angular distribution is shown in Figure 4. The conductor is a flat copper wire 
doubly wrapped with glass and has a cross section of 1.8 x 2.65 mm. According 
to actual measurements,‘ the uniformity of the bipolar field is better than 

1 percent when R < 15 mm. The total inductance of the deflector coil is 257 


BH . 


If the uniform bipolar magnetic field B generated by the deflector coil rises 
linearly with time, then the following relation exists when the e-beam is 
deflected by a small angle 6 (Figure 5): 





200/E* + 2 EE, (8) 
31(L + 1/2) 





where E is the kinetic energy of the e-beam (MeV), E, = 0.511 MeV, L is the 
distance from the exit of the deflector coil to the fluorescent target (m), Y 
is the vertical coordinate of the scanned e-beam on the fluorescent target 
(m), and 1 is the effective length of the deflector coil. Actual measurements 
showed that 1 = 0.2032 m and L = 1.8783 M. In order to have the image of a 35 
MeV e-beam fall at Y = 25 mm, the deflecting magnetic field should be: 
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Figure 6. Magnet Excitation Curve of 
Figure 5. Relation Between the Sweep Coil 


Drive Current and Deflecting 
Power of the Sweep Coil 


Figure 6 is the measured magnetic excitation curve of the deflector coil. It 
is obvious that the excitation current required is less than 30 A. The design 
specification of the magnetic excitation pulsed power supply is that the pulse 
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current amplitude is +30 A. The leading edge rises linearly and the rise time 
is 2~5 ws, adjustable at different scales. This power supply is being 
developed. The maximum inverse peak voltage of the deflector coil is: 


pat - 257 x 10-6 x —82 


which meets the voltage breakdown requirement of the amplifier of the power 
supply. 


Because the vertical "object slit" is designed to be 1/6 mm in full width, the 
effective working width of the fluorescent target in the vertical directionis 
50 mm. Hence, within a macropulse, it is possible to obtain 30 points in the 
time-dependent energy spectrum for diagnostics. For a macropulse which is 2~5 
us wide, the camera is required to have a resolution of better than 60~150 ns. 


3. Fluorescent Target, Video Camera, and Computer Graphics Processing System 


An imported French-made Cr-activated Al,0, fluorescent target (AF995R) with an 
effective working area of 100 x 50 mm was used. The TV video camera used is a 
Model SGB-11N, which has 350 vertical lines (<0.2 mm) and 550 horizontal lines 
(<O.2 mm) and a time resolution of approximately 8 ns. The camera is equipped 
with a microchannel-plate light intensifier and the amplification factor is 
10,000. Further, the applied bias voltage between microchannels can be used 
for high-speed shutter control. The shortest gate control interval is several 
dozen nanoseconds. Therefore, it can be used to take the energy spectra of the 
e-beam at different intervals within the same macropulse in order to study 
imortant mechanisms such as the buildup process of the FEL oscillation and the 
emission efficiency of a micropulse. Figure 7 shows a block diagram of the 
image processing system. 
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ing figures, Figure 7. Block Diagram of Image Processing 


tables, and 
data can be provided on-line by the energy spectrometer by using the image 
processing software and FEL cesearch software. 


a. Macropulse time-dependent energy spectra (current vs. energy, current 
vs. time, energy dispersion vs. time, instantaneous current vs.time, 


etc.). 


b. FEL emission efficiency. 
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[Text] Abstract 


An analytical study is done on mode locking in an FEL (free electron laser) 
and the mode-locking conditions are obtained. It is derived that due to mode 
coupling the number of active modes in the cavity is reduced by a factor of m, 
the low gain factor is increased by a factor of m, and the exponential growth 
factor goes up by a factor of m,’/?; m is the ratio of the micropulse 
repetition frequency Av, to the frequency interval of modes in the cavity Ay, 


i.e.,m = Av,/Av,. 
I. Introduction 


In an RF-Linac-driven IR short-pulse FEL, when the micropulse length lk is 
comparable to the total slip between the wiggler and the light beam A, i.e., 
lh = A, the propagation of the e-beam, i.e., the nonsteady-state term, must be 
taken into account in the equation. The mode locking under such conditions has 
been investigated by numerical modeling in reference [1]. However, when the 
length of the micropulse, Iy,, is still much larger than the total slip between 
the wiggler and the light beam A, i.e., ly >> A (such as ry~ 10 ps, A,~l po, 
Lp = ¥zTp = 0.3 cm, A = NA, = 100 x 10°“cm = 0.01 cm), the slip effect between 
the e-beam and the light beam can be neglected. It can be treated as a steady 
state. The mode-locking problem in this case can be studied analytically. This 
is an attempt to analytically study the mode-locking conditions and mode- 
coupling gain of an IR short-—pulse FEL. 
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II. Gain of Coupled Mode IR Short-Pulse FEL 
1. E-Beam Pulse Profile 


Let us assume that the e-beam pulse follows a Guassian profile, i.e., 





n(to) = Nof (to) (1) 
— 1 Vglp . -v2e2/20%, 
f£(t,) = x 4, e (2) 


where no, vz, 0, and Tp are the density, longitudinal velocity, mean square root 
of the spatial length of the pulse, and time width of the pulse, respectively. 


2. Optical Mode Evolution Equation 


If the e-beam pulse length is still greater than the overall slip between the 
e-beam passing through the wiggler and the light beam, the problem can still 
be treated as a steady-state case. According to reference [2], the equation 
for the FEL longitudinal mode which takes the pulse structure of the e-beam 
into consideration is as follows: 


With a linearly polarized wiggler, we have 








da,(Z) _ wa, -idk,z 1 2 iw Co 3 
re Co 7 ” 
a, = -ai se! (4) 
V2 mc*k,, 
a: (5) 
V2 me*k, 
2 
wi, . 4ne*n, (6) 
m 
Fg = Jo (€) - Ji (&) (7) 
2 
E = —_—* (8) 
2 (1+a,) 
1 
AKn = Ky ~ Kal ioy - 1) (9) 
On | 2kYo (10) 
Cc yu? 
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where E,, w,, and k, are the complex optical field, frequency, and wave number 
of mode n, respectively; e, m, and ng are the charge, wass, and beam density 
of the electron, respectively; B, and k, are the amplitude and wave number of 


the wiggler field; and J, and J, are the zeroth- and first-order Bessel 
functions of the first kind. 


In an FEL, due to the structure of the e-beam, mode coupling can occur in the 
linear region. Hence, 








_ _ io, -10 gt, ei - 1 ,, -10 ,t, f# 1 
F,(Z,€,) = exp{-——* ye [Lnan} 1 = ye [ bnandz (11) 
2 /a @® “ 
L - a 7 rs) m {Ak ez Ni 
. wed | 7. F, ( @, ©) kne dz (12) 


By substituting equations (11) and (12) into (3), we have 


T 

















da, (2) - Wray (6). -iak,z 1 ? dept (ty) eta 
dz 4cy,®, T,J-+ 
2.2 Tp 
F 2k.)? _; 1 Co 
. i W pay (6) ( : af e ‘ako  f 7,dt,f(t,)xe ve ws re (o. 
4cy, , T,/-> . y Fy ( - 6) 
x [‘ae’ {* dz"eith a (2) (13) 
0 0 
Let 2 
e 1 2 1@ ,Co 
f, = ae dt,f(t,)e (14) 
T» 
1 2 i(w,-@,) Cy 
f, = Tl dt,f(t,)e (15) 
By substituting equation (2) into (14) and (15), we have 
f, o in! s/2V, (16) 
f, . eo (MnO _)20%,/2v, (17) 
For an IR short—pulse FEL, 
=f « 2 o 10*(1/cm) 
v; c 
L 
0, = —2 ~ 10°'cm 
—— tune)? @ 10° 
z 
. f, = 0 (18) 
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By substituting equations (17) and (18) into (13), we have 


_ (@,-@,) 29% 








da, (2) _ , @ranFy(€) (2k) 2e7 ‘i re ™ 
~ Oz 4cy,@ ,h? = 
x —#F (—*) ['dz! [* dz"e"**'a,(z") (19) 
n n 


Equation (19) is the mode-evolution equation of the IR short-pulse FEL. As we 
can see, if the following inequality is valid: 


7 2g? 
(@ ,-@,) “Oo, 


Vv 3 


<1 (20) 





then it is possible to couple several modes. From equation (20), the pulse 
length of the e-beam is dependent on the number of coupled modes. This is 
because the farther apart wp, from w, is, or the larger the number of coupled 
modes, then the shorter o, is required to be to keep equation (20) valid. 


3. Number of Coupled Modes 


According to reference [1], the number of modes that can be simultaneously 
activated in the resonant cavity is 


Av 
~%° By. (21) 
where Avg is the gain bandwidth. It is at least determined by the uniform-gain 
bandwidth and Avg = v,/2Ny, where vg is the center frequency of the light beam, 
N, is the number of periods of the wiggler field, and Ay, is the frequency 
interval between modes in the resonant cavity; Av, =c/2l,, where Lk. is the 
length of the resonant cavity which is at least equal to lL, and L, is the 
length of the wiggler. 








. M&> og He tty (22) 
& 


Therefore, thousands of modes may be simultaneously activated in the cavity. 


In reality, the number of activated modes decreases with growing coherent 
length. Nevertheless, this kind of growth is limited by the pulse length. The 
shorter the pulse length, the larger the number of activated modes. 


As described in reference [1], an interferometer inside the cavity (e.g., a 
slightly oblique etalon) might be used to turn independent pulses into 
coherent light pulses. The number of light pulses that might be coupled is 
determined by the number of independent coherent pulses. The number of 
independent pulses is further determined by the repetition frequency of the 

e-beam pulse. Since we adjust the length of the etalon with respect to the 
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repetition frequency of the e-beam pulses so that the pulses in the etalon 
overlap with the new ones, the number of modes to be coupled can be determined 
by the following equation: 





m,. = z (23) 


where Av, is the repetition frequency of the electron micropulses. It is given 
in reference [1] that L = 6m, Av, = 1 GHz and m = 40. 


From the number of coherent pulses, m = 40, the maximum difference between the 
frequency of the coupled mode and the given laser frequency w, can be 
calculated. By substituting this difference and the given electron pulse 
length of o, = 0.1 cm into equation (20), it is found that the equation is 
valid. Hence, m = 40 is the number of modes that can be coupled. 


4. Low-Gain Situation 


Let a const (24) 


lja,|e’** 


a, = |a,|e’** = const (25) 


By substituting equations (24) and (25) into (19), we have 
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4c*y, 
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(1) Mode-Locking Conditions 
~w .) 20? 
(@ ,-@ ,) *O, 3 
2v4, 
{= 9, (30-1~4) 
la,| ~ 4, 
6, ~ 8,(0, - w,) 


(n= 1, 2, -:+, m; m=1, 2, ---, m) 





Equation (30) indicates that when the phases (and amplitudes) of m short-pulse 
modes of almost the same frequencies are identical, they are coupled into a 
single optical mode. 


(2) Gain of Coupled Mode 


Using equation (30) in (26), we have 








G. = mG, (31) 
‘.. gd , sin*é 
G, - - 2G, (>)” a8, | Q2, ) (32) 
22.2 
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(3) Power of Coupled Mode 
5. High-Gain Situation 
a, = a,e7* at*et (35) 
Let 
a, = a, eV orte (36) 


By substituting equations (35) and (36) into (19), we get 
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(2) Gain of Coupled Mode 
By applying equations (37) and (38) to (39), we get 


2 (m.G,) (Ak,+® ,) a, 
a, = (40) 
[ (Ak,+® = *+a7] . 





- 2»? 
= - (mG) AK, ® )?-a4) (41) 
[Ak,+® _)2+a4)? 





From equations (40 and (41), we can obtain the following solutions: 








- 3 G y1/3 le Ak, 2 4? 

a, o (m.G,)*/3[ (Sime, ] (42) 
- i 1/34 - Ak, 

®. > (m.G,)*/? [1-4 ; (m.G,) 7 (43) 


(3) Power of Coupled Mode 


P. = P,exp{y3 (mG) /L, [1 - ( a __)2)} (44) 


3 (m.G,) 1/3 





III. Conclusions 


The gain of the coupled mode of an IR shoi -pulse FEL was studied and mode- 
locking conditions (equations (30) and (39)) were obtained. Due to mode 
coupling, the number of active modes in the cavity is reduced by a factor of 
m. and the low gain is increased by a factor of m (equation (31)). The 
exponential growth rate will rise by a factor of m,’/*? (equation (42)); m can 
be determined by equations (20) and (23). 
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[Text] Abstract 


The distribution of coating thickness across a 0.8-m-diameter mirror is 
discussed under the geometric constraints of the coating chamber. The effect 
of nonuniformity of the metal-enhanced high-reflective-index coating on the 
polar’ zation of the light beam is presented. Phase distortions associated with 
tne 1.06-ym laser and the 0.5893-ym beacon light are calculated. By means of 
design modification, phase changes of the same sign are obtained to allow the 
beacon light to achieve its expected objective. 


I. Introduction 


The diameter of the mirror on the emitting end of a high power laser is D. 
Consider the light beam following a Gaussian distribution; the angle of 
divergence at the waist is § = \/D (where A is the wavelength of the laser). 
Furthermore, considering other nonlinear optical effects such as the thermal 
blooming caused by the transport of high power laser through the atmosphere 
and excited Raman scattering, the mirror diameter D ought to be large. In 
practice, the main mirror of a high power laser is usually 5~20 m in diameter. 
A large diameter is difficult and costly to make. Moreover, the preparation of 
its coating, especially the uniformity of the coating thickness, is extremely 
critical and difficult. The nonuniformity causes polarization aberration and 
wave-front distortion which affect the far-field diffraction characteristics.? 


Another important issue of a high power laser system is that there are many 
factors causing wave-front errors in the optical path. The energy of the laser 
beam is dispersed, which needs to be compensated using the adaptive optical 
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technique. However, in order to make this kind of correction possible, a 
suitable light source, usually called the beacon light, is required as the 
information carrier for the wave-front error. 


A phase shift is produced as a light beam is reflected by a coated mirror. 
Different phase shifts are generated by light of different wavelengths. 
However, due to nonuniformity of the coating thickness, phase distortions are 
generated. At different wavelengths, the phase distortions are also different. 
It is even possible to have phase distortions of opposite signs. Hence, the 
opposite effect will be produced when it is compensated by the beacon light. 
Thus, the expected result cannot be achieved. 


Based on existing conditions, the coating thickness distribution on a 0.8-m- 
diameter mnear-spherical mirror is calculated by taking the geometric 
constraints of the coating chamber into consideration. Furthermore, this 
distribution is optimized. The phase distortions produced by this type of 
metal—enhanced high-reflectance coating series at the laser wavelength, i.e., 
1.06 wm, and the beacon light wavelength, i.e., 0.5893 wm, are analyzed. By 
modifying the design of the coating, the signs of phase distortions at these 
two wavelengths are consistent. 


II. Effect of Physical Constraints in Coating Machine on Coating Uniformity 


Rotating 
; : sos _ . primary 
There are two types of uniformities concern Mirror Ae 
ing the coating of a large-diameter mirror, diameter D, \ 
—_—— ound 





i.e., orientation uniformity and _ radial 
uniformity. Orientation uniformity can be 
easily compensated by spinning the mirror at 
a high speed via an evaporating ring such as 
the one shown in Figure 1. Usually, for a 
large-diameter mirror, from a safety stand- 
point, a low spin rate is desired; hence, it 
is necessary to use a ring-shaped vapor Source-to- 
source to achieve excellent orientation oe peor 
uniformity. However, radial distribution of 
coating thickness is more difficult to con- 
trol. Radial distribution of coating thick- 
ness depends upon three parameters’: 1) the F 
value of the mirror and the nonspherical 
shape of the surface, 2) the distance between 
the vapor source and the substrate H, and 3) 
the source offset D from the vapor source to Vapor source 
the center of rotation of the mirror D (off-— Figure 1. Ring Source Coating 
set). With a spherical rotation fixture, the Geometry Using Series of Vapor 
coating thickness at an arbitrary point P on Sources Arranged to Simulate a 
the substrate is Ring Vapor Source of Diameter 
2D Placed Underneath a Rotat- 
ing Primary Mirror 
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[h2+ (D+p) ?) 3/2. 
where 


h = H —- R (1 — cosA) 


sina = £ 
R 


. mb , [(h*+D*+p*) (hcosA+psinA) -2pD*sinA) (1) 
[h?+(D-p) ?) 3/2 


(2) 


(3) 


In the above equations, m is the total weight of the evaporated material, uw is 
its density, R is the radius of the spherical mirror, and p is the radial 
radius from point P to the axis of rotation. The relative thickness is the 
ratio of t,, to t,,, which is the thickness in equation (1) when p = 0. 


Figure 2 shows the optimal condition for an 
0.8-m—-diameter F/2.0 mirror to achieve coating 
uniformity. From this figure, one can see that 
in order to achieve better than 0.5 percent in 
uniformity, the distance between the vapor 
source and the substrate needs to be 0.9 m and 
the offset is 0.63 m. The coating machine 
should have a height of 1.3 m and a diameter of 
1.5 m. However, our equipment does not meet 
these requirements. Due to cost considerations, 
modifications will have to be made for the 
existing equipment. Therefore, it is difficult 
to achieve the ideal coating uniformity. Figure 
3 shows the effect of offset D on the coating 
uniformity at a known height H. From the fig- 
ure, when H is 0.65 m, the offset should be 
0.48~0.5 m to achieve a uniformity of +1.5 per- 
cent. Nevertheless, the diameter of our exist- 
ing coating machine is 1.2 m and the maximum 
offset one can attain is 0.46 m. Figure 4 shows 
the coating thickness as a function of height 
with a bias of 0.45 m. This figure shows that a 
+2 percent uniformity may be obtained at a 
height of 0.56~0.64 m. 
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Figure 2. Optimized Coating 
Thickness Uniformity for 
0.8-m—Diameter F/2.0 
Concave Mirror 


In the above discussion, shields are not used. In all uniformity calculations, 
all vapor sources are small-area sources and obey the cosine distribution. The 


mirror is considered to be spherical. 


III. Initial Design of Coating Series 


The main laser wavelength is selected at near 1 wm, e.g., 1.06 wm. The beacon 
light might either be scattered laser light reflected by the target or a laser 
light of another wavelength on a relay satellite. We are using the artificial 
beacon light produced by the resonant scattering of sodium atoms in the alkali 
metal layer formed by falling meteorites 80~100 km high in the atmosphere in 
response to an 0.5893 wm laser as proposed by French astronomer Labyrie. It is 
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formity for 0.8-m—Diameter F/2.0 Figure 4. Coating Thickness Nonuni- 

Concave Mirror for Various Offset formity for 0.8-m-Diameter F/2.0 

Distance and for a Fixed Source- Concave Mirror for Various Height 

to-Substrate Distance of 0.65 m and a Fixed Offset Distance of 
0.45 m 


assumed that a reflectance of greater than 99 percent is required at 1.06 pm 
and that at 0.5893 wm it needs to be higher than 95 percent. 


The initial design of the coating is 
G/Cr/Ag/A1,0,/(LH)*/Air 


where L and H are the one-fourth-wavelength thicknesses of SiO, and Ta,0. at 
the laser wavelength, respectively. Cr and Al,0, are bonding layers. G is the 
substrate and the geometric thickness of the Ag layer is 70 nm. 


The design with a full or half one-fourth-wavelength dielectric film on a 
metal layer is selected because it is easy to monitor. Initially, the (LH)? 
structure was chosen because it has wide bands of reflectance near the two 
wavelengths of interest. Ag was selected because it has a high reflectance 
band ranging from visible to infrared. Based on this design, the maximum 
incident angle of the main mirror is 15°. The following is a discussion of the 
effect of coating nonuniformity on the polarization characteristics (elliptic- 
ity and phase distortion) at the two wavelengths of interest. 


IV. Effect of Thickness Nonuniformity on Ellipticity of the Initial Coating 
Design 


When a linearly polarized light shines on the main mirror at an oblique angle, 
the reflected light is elliptically polarized. In other words, there is a 
phase difference between the vector component P which is parallel to the 
incident plane and the vector S which is perpendicular to the incident plane. 
Furthermore, the perpendicular components have different amplitudes. Figure 5 
shows the changes in reflectivity and phase as a function of thickness 
nonuniformity at the main laser wavelength and beacon light wavelength. At 
1.06 wm, the phase difference Ap-As < 1°. Even with a 4 percent thickness 
nonuniformity, the difference in reflectivity is less than 0.02 percent. 
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At 0.5893 wm, the phase difference is less 
than 3.5° and the reflectivity difference 
is less than 0.2 percent. In conclusion, at 
an incident angle of 15°, the reflected 
light is elliptically polarized in a very 
flat manner. A 4-percent thickness nonuni- 
formity has very little impact on its 
ellipticity. 


V. Effect of Coating Nonuniformity on Phase 
Distortion of Reflected Waves at Different 
Wavelengths 


The characteristics of nonuniformity of 
thin coating include different optical 
thicknesses and phase shifts (phase distor- 
tions). Phase shift will cause image dis- 
tortion. In particular, it significantly 
affects the far-field diffraction charac- 
teristics of the system. Figure 6 shows the 
phase distortions of the reflected waves at 
the two wavelengths of interest as a func- 
tion of coating thickness nonuniformity for 
the initial design. As one can see, at 1.06 
um, a +4 percent nonuniformity only results 
in an -11° phase shift or a wave-front 
distortion of -\/32.7. A 1-percent nonuni- 
formity causes a phase shift of -2.4° 

or a wave-front distortion of 4/150. At 
0.5893 pm, the phase distortion is 
larger. A +4-percent nonuniformity 
results in a phase shift of 140.28° or 

a wave-front distortion of 4’/2.57. A l- 
percent nonuniformity causes a phase 
shift of -8.5° or a wave-front distor- 
tion of A/42.3. 
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Figure 7 shows the radial distribution 
of phase distortion at the wavelengths 
of interest for an 0O.8-m-—diameter, 
F/2.0 mirror. The calculation is based 
on an assumption that the large mirror 
surface is close to spherical. The 
distance between the vapor source and 
the substrate is 0.55 m and the offset 








=F =7 =3 7 2] 3 4 
r— 10 





Figure 6. Phase Distortion in Degrees 
for the Initial Selected Design at 


is 0.45 m the Two Wavelengths as a Function of 
. . Thickness Nonuniformity for 0.8-m- 


Opposite phase distortion signs cause 
other problems as well. An optically 


Distance F/2.0 Mirror 


erroneous image at the wavelength may be partially corrected by means of self- 
adaptive optical compensation. If the wave-front distortion of another 
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wavelength has the opposite sign, the 
correction at one wavelength will 
cause the further deterioration of the 
wave front at the other wavelength. In 
order to simultaneously correct the 
wave fronts at two wavelengths, phase 
distortions at these two wavelengths 
ought to have the same sign. To the - 
extent possible, the phase distortions 
of these two wavelengths must be con- 
sistent. To this end, a modified de- 
sign was adopted. 


VI. Modified Design and Its Effect on Fi 
Polarization and Phase Distortion S 
t 
To have the same phase distortion f 
signs, numerous simple designs were 0 
investigated. Based on _ theoretical 
calculation, the modified coating was des 


G/A1,0,/Ag/A1,0,/ (LH) 22 /air 


H and L are the one—fourth-wavelength Ta,0. 
and SiO, correction coatings for 1.17 pm, 
respectively, and Al,0, is the bonding 
layer. The effect of thickness nonuniform— 
ity on the polarization characteristics of 
this modified design is shown in Figure 8. 
It was found that the reflectivity of this 
design varies significantly with thickness 
nonuniformity at an 0.5893-ym wavelength. 
The two designs have comparable reflectiv— 
ities at 1.06 wm. Compared to the initial 
design, the reflectivity of the modified 
design is slightly higher than a +4-per- 
cent nonuniformity. At 0.5893 pm, the 
reflectivity of the modified design is 
lower than that of the initial design. 
Figure 9 shows the phase distortions 
produced by the modified design at two 
wavelengths. It shows that a 4-percent 
nonuniformity results in a -23° phase 
shift, or a A/15.6 wave-front distortion, 
at 1.06 pm. A 1-percent nonuniformity 
causes a phase shift of -4.5° or a 4/80 
wave-front distortion. At 0.5893 pm, the 
distortion is larger: a 4-percent nonuni- 
formity produces a -55° phase shift, or a 
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Pe A= 0. 6893 um 
30 l—percent nonuniformity results in a 
-14.5° phase shift, or a ’A/24.8 wave- 
P— 40 front distortion. The wave-front distor- 
tions corresponding to different coating 
Figure 9. Phase Distortion in thickness nonuniformities have identical 
Degrees for the Modified Design signs at these two wavelengths. 
at the Two Wavelengths as a 
Function of Thickness Nonuni- Figure 10 shows the radial distribution 
formity for 0.8-m-Diameter of phase distortion for the modified 
F/2.0 Mirror design at the two wavelengths for the 


0.8—-m F/2.0 mirror. 
VII. Conclusions 


The analysis described above indicates that when using the modified design 
(adding one-fourth- or one-eighth-wavelength dielectric coatings on top of 
metallic silver), coating thickness nonuniformity could produce the same phase 
distortion signs at two wavelengths. At an incident angle of 15°, the modified 
design does not have much impact on its polarization characteristics. The 
reflectivity in the wavelength range of interest does not decrease. This is a 
suitable design for an 0.8-m-diameter mirror. For ease of monitoring, all 
dielectric coating was controlled by a single wavelength to minimize any 
variation caused by dispersion due to the use of different wavelengths for 
controlling different layers. 
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[Text] Abstract —— 


On the basis of work done on vacuum synthetic-aperture beam control, the 
effect of phase control errors (including tilt, piston, and focus errors) on 
coherent beam combination is analyzed and discussed. A corresponding physical 
model is presented. In addition, numerical computation was done and is illus- 
trated in the form of three-dimensional plots. 


I. Introduction 


The use of a synthetic aperture as a coherent high-power laser emission system 
has attracted a great deal of interest in recent years.’ By controlling the 
wave-front phase of the light beams from each subaperture of the synthetic 
aperture system, it is possible to achieve coherent combination of the light 
beams in the far field. Nevertheless, phase-control errors usually occur with 
wave-front phase control. There are several factors affecting the control of 
phase which leads to a decline in performance of the synthetic aperture system 
and lowers the far-field peak of the combined beam. We have studied the beam- 
control concept under ideal conditions (i.e., without various factors 
affecting the precise phase control of the beam). A physical model for the 
alignment, transition and coherent and noncoherent combination of the beams 
was presented. Furthermore, theoretical analysis and numerical computation 
were performed.*:? On the basis of our previous work, this paper discusses 
factors affecting the precision of phase control (e.g., piston error, tilt 
error, and focus error) and their effect on coherent beam combination. 


46 











A physical model is presented and numerical computation is carried out; 
results are illustrated in the form of three-dimensional figures. 


II. Theoretical Analysis and Numerical Computation 


To facilitate the theoretical analysis, let us assume that the laser emits a 
planar beam. 


Up(Xy, Yi) = Bpy (X,, yy) (1) 


where B is a constant which is there for use in the analysis to follow. Let 


B = 1 at this point; p,(x,, y,;) is the pupil of the light emission system which 
is defined as follows: 


(2) 


1 xi+yi<(d/2)? 
P,(X,, Y¥,) = 


0 (xf+y?)>(d/2)? 
where d is the diameter of the pupil. 


Let us assume that the synthetic aperture is composed of N = 1 x m apertures 


and the center coordinate of each aperture is x,(i,j), ys(i,j),and i-1,2,---,1l 
and j=1,2,---,m. The coherent combined far field is?: 
1 oo (x-x,) /Az 
U(x, y) = A(x,y,z) Ff > P)(7,,5)}] 
»» sy I" y-y,) faz 
(x-x,) /Az 


A(x, y,z) F{p, (r,)} 
trie My-y,) /Az 








1 | X-X, | 7 - 
p> exp{-2nj7[ le x, (i,j) + iz y,(i,j))} 


where A(x,y,z) = exp (jkz) - exp (( DE) (x@sy*)1/( jaz), F({ )} represents the 


Fourier transform, k = 2n/A where A is the wavelength, and (x;,y,) is the 
coordinate of the target. Let xy = yy = O (i.e., all the beams from the 
subapertures converge on the axial axis), then we have 


x/Az 
U(x,y) = A(x, y,z) F{p,(1,)}| 
y/az (3) 


-2nj(% i 4 Sy (ij 
? exp { 2nj(>—x,(i,J) + pe ¥elt FI} 


1 
asl J 
Obviously, from equation (3), the peak light intensity of the coherently 
combined beam is: 


Iy = |U(0,0) |? = gage Ftp, (14) ) 2 * (1xm)* = N* + Iiy (4) 
Zz 


where Ij is the peak intensity of each individual subaperture system at 
(0,0,z). 


4/ 








In the following analysis and computation, let N = 16, d = lm, A = 1 pm, z = 
10°km. 


Figure 1 shows the intensity distribution of the coherent combination of beams 
by the synthetic aperture. (The numbers shown are relative values.) Figure 2 
shows the energy distribution curve in the far field; Rp represents the Airy 
disc radius: E(AIRY DISC)/E(TOTAL) = 0.868299E+00 indicates that the energy in 
the first airy disc is 86.8299 percent of the total energy. Rp = 0.61 x A x 
z/b, which is in units of 0.61 x A z/b where A is the wavelength, z is the 
distance from the mirror to the target, and b is the effective pupil radius of 
the phase-control array. 
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Figure 1. The Intensity Distribution 
of Coherent Beam Combination E (AIRY DISC)/E( TOTAL) = 868 % 
Figure 2. The Regional Energy 
Distribution in the Far Field. 
1. Effect of Piston Error in Phase E = 86.8% 
Adjustment‘ 


When performing wave-front phase adjustment and control of a divergent beam 
from a subaperture, a phase error in the form of exp (jkgjj), where 9j; is 
independent of x, and y,, is the piston error. Hence, the syntheticaliy 
combined far-field distribution becomes: 


x/Az 
U,(X,¥) = A(x, y,2) Fip,(1r,)}| 
4 . y/Az 
(5) 
] m 
“SY expl-2nj(~ x, (i, 5) + Soy,(i,5))] + exp (jk ,,) 
I*l jel Az ° Az 
Consequently, the coherently synthesized peak intensity is 
Lup Y(0,0) | 
. lo | 1 ‘ (6) 
-——| F{p,(r,)} exp (j7k®@,,) 
12> ie » ss P\J ® 
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neral., 5 — (jk yl N?, 
»» > exp (jk@, < 


I wp < N°: Ty 


that 


1is means the piston error would cause the peak intensity to 


fall. 
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Figure The Intensity Distribution With the Piston ag = ) 
; pede E (AIRY DI$€)/E( TOTAL} 
(O,=%,,"=%,=9,=0.51. other g = ) ) 
oa: . er @,,20) Figure 4. The Regional Energy Distribution 
in the Far Field With the Piston 
Errors 
(9,79 =, =@, =90.5A, other 
Figures ind 4 show the far-field light ? ) 
eure’ ele aane- 4s o =0), E=86.5% 
intensity distribution and energy distri- 
but Lor ré spectively, with Wi. = ~12 ™ Yo 
Y2: \ and @jj = 0. These figures show that piston error will lead to 
decline of the peak intensity of the combined beam. Howeve1 it has a 
relatively small effect on the energy in the first Airy dis: 
2. Effect of Tilt Error‘ 
is a phase error term in the form of exp({|jkKicjj,x(x-x,(i,j)) + 
VC y¥—-Y¥e(i, j which is caused by the tilt and vibration of the subaperture 
his case, the optical field of the combined beam can be expressed as 
] m x/Az 
x, x,y, z,)¥ yo F{p, (r,,;) expljkc,,,(%,-X,) +7kKc;4,(¥,-V,~) ]} 
fF, v/Az 
os x/hz 
A(x, Y,2, » F(p,(r,) exp (JkE 45%, 4C;5¥,) 1} ] 
1 y/Az 
ext A ) x x + y y ) ) 
- Az * Az *® 
loom (X-CyyqZ) /AZ 
4 xy 2S F{p,(r,)} . 
fe? je] (x=-C jy= Az 
~ ; x y 
ex] nT X,(1,Jj)+ ' 








It is obvious that the peak intensity becomes 





lu = U.(0,0) |* 
1 om -¢ /d ] ~ (8) 
l = tie | 1 0 
. =e F{p, (r,)} fe >» F(p,(r,)}|. |? = Nz, 


Figures 5 and 6 are the intensity and energy distribution when cjjx = cjjy = 0.5 
x 10 °rad. Based on these figures, we can see that tilt error not only causes 
the far-field peak light intensity to fall but also results in a reduction of 


energy in the first Airy disc. The larger the tilt error, the more apparent 
this effect. 
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Tilt Errors. (¢ =¢ -=0 4S» 10 ‘rad ) 


Figure 6. The Regional Energy Distribution in 
the Far Field With the Tilt Errors 
3. Effect of Focus Error (c =e, = 05x 10°*rad) E = 85.6% 


When focusing the beam, if it is not 


accurately focused, a focus error term such as the following will be intro- 
duced?: 


exp (jkejj[xy-xs(i,j))? + (yi-ys(i,j))?7)/2) 


Hence, the far-field combined beam becomes: 


] m xX Zz 
Up(x, Vv) = A(x, y,z) > F{p, (r,,,) exp (jke,,((x,-x,)*+ (y,-y¥,)*) /2)} v/he 
i=l je] é 
—. x/Az 
= A(x, y,z) } YO Flp,(r,) exp (jke;;r;/2)} 
o—< j=7 : y Z 
- X ,. 
SOAP USB le Ags ses!” po ¥ekt J ) } 
4 a 
Let € 4, . , the 
) : ” 
x \ a xpiyAA re 22) 
viAaz 
aad od, ~ 4 (9X 
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Hence, the peak intensity of the combined beam is 








‘ > ] : = . ‘ 0 
| = |\U,(0, O)}° = j252 > F{p, (r,) exp (jrA,,ry/Az) } 0 
Z°| ied 31 ; 
(10) 
1 . 0 | 
_ = N? 
122? (2. »» Fip,(r,)} 0 N* Iie 








Figure / is the intensity distribution of the coherently combined beam at z = 
5 x 10°m without any focus error (Ajj = 0). Figure 8 shows the intensity 
distribution at z = 5 x 10°m with a focus error, i.e., Aj; = Ajz = Az, = Ay = 
A3, = 432 = Agi = Ag2 = O.1 and Aj3 = Ay, = O23 = Azy = Az3 = Azy = Aug = Any 

-Q0.1; the physical significance is that some of the subapertures are focused 
in front of the common target and others are behind it. It is apparent that 
the peak intensity drops considerably in this case. Moreover, the larger the 
focus error, the more severe the impact on the beam combination. 





Figure 7. The Intensity Distribution of Coherent Figure 8. The Intensity Distribution With the Focus 
Beam Combination Without Focus Errors Errors(A =A =A =A_-A #4. #4 
(3 =0).2=5x 10° m. 4.70.1, 4,-4,=4,-4, =a =A 2d. 


=A.=-O01). 2=5«10'm 


III. Conclusions 


A detailed analysis has been done on factors affecting the accurate control of 
the phase of beam combination by means of a synthetic aperture; furthermore, 


numerical computation has beer ne. The results show that the peak intensity 
of the combined beam drops c jerably because of failure of accurate phase 
control. However, it shovid vinted out that the factors considered in the 
analysis were assumed to be atic, instead of random and dynamic variables. 
As for laser transport in the atmosphere, it is more meanir to investigate 


the effects of random and dynamic variables, such as cuoulence, on the 
accuracy of phase control. We plan to work in this area in order to establish 
a realistic model to study the effect of turbulence on phase-control accuracy. 
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[Text] Abstract 


The emission characteristics of velvet cathodes have been experimentally 
studied under long (50~100 ns) and short (3 ns) pulse conditions with 
corresponding cathode densities of 400 A/cm* and 7 kA/cm’, respectively. 


The current-voltage curve was measured. The cathode emission pattern was 
obtained using the Cerenkov effect. Scanning electron micrographs of the 
velvet cathode were obtained and the experimental structure of the velvet 
cathode was measured. The emission mechanism of the velvet cathode is 
discussed. 


I. Introduction 


An intense relativistic electron beam (IREB) has a very high current density 
and a very low temperature. It is a prerequisite for a number of devices such 
as electron-beam—pumped excimer lasers and FELs. Different applications have 
different characteristic requirements for the e-beam. An e—beam—pumped excimer 
laser requires a beam with very good uniformity and the plasma between the 
anode and cathode should have a low closure speed. A Raman FEL, on the other 
hand, requires a high brightness and low emissivity. 


These requirements can be met by using different diodes and cathode materials. 
Recently, it was proven that velvet is an excellent cold-cathode material. Its 


major advantages are as follows: First, its start-up time is very short,* i.e., 
its impedance can drop from infinity to its characteristic value instanta- 
neously. Second, the anode-cathode plasma closure speed is very slow and 





plasma expansion is almost unnoticeable. Thus, the current profile of a diode 
with a velvet cathode is almost an ideal square wave. Third, compared to other 
cold-cathode materials, the velvet cathode is several orders of magnitude 
higher in brightness. 


However, the emission mechanism of the velvet cathode is still not well 
understood. To this end, some experimental work is done in order to study the 
emission mechanism involved. 


II. Experimental Apparatus and Results 


In this experimental work, three different IREB accelerators are used; their 
pulse widths are 120 ns, 50 ns, and 3 ns. The 3-ns device is called the 
Febetron 706. They are essentially composed of three parts: 1) An energy 
storage unit, which is a Marx generator for the 120-ns and 3-ns apparatus. The 
50-ns apparatus uses a high-voltage transformer as its energy storage unit. 
2) A pulse formation line, which is a water Blumlein line for the longer-pulse 
apparatus and an oil-based Blumlein line for the Febetron 706. 3) A diode 
which has a rectangular e-beam cross section for the longer-pulse apparatus 
and a circular e-beam cross section for the short-pulse device. The character- 
istics of these IREB accelerators are tabulated in Table 1. 


Table 1. Specifications of the IREB Accelerator Systems 




















Model IREB source (1) IREB source Febetron /06 
Parameters (I1) 
Output voltage (kV) 300 250 600 
Output current (kA) 30 25 6 
Output impedance 10 10 100 
(Q) 
Pulse width (ns) 120 50 3 
The shape and area Rectangular Rectangular Circular 
of cathode 40 x 4O cm 20 x 3 cm ~9 mm 























The results are as follows. 
1. Measuring Cathode Emission Pattern Using Cerenkov Effect 


Real-time measurement of the e-beam current was measured by means of Cerenkov 
radiation using an optical fiber or a quartz plate.* This measurement is real- 
time and is not susceptible to interference. 


This experiment was carried out using IREB accelerator (II). The experimental] 
set-up is shown in Figure 1. Figure 2 [photograph not reproduced] depicts the 
emission pattern of the velvet cathode. It clearly shows that the velvet 
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Figure 1. Schematic of Experimental 
Arrangement 


cathode emission is more uniform. Its 
emission pattern consists of a number 
of bright and dark spots. The sizes 
and locations of these points seem to 
follow a certain pattern. 


The structure of the velvet cathode 
was observed using a microscope and 
is illustrated in Figure 3. The vel- 
vet fibers are essentially standing 
straight, and are approximately 20 pm 
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Figure 3. Practical Construction of 
Velvet Cathode (Unit: mm) 
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in diameter. There are a number of gaps at the back side of the velvet cathode 
and the sizes and locations of these gaps also follow a certain pattern. We 


measured the sizes of all the bright and dark spots. 


consistent with the sizes of the velvet cathodes. 


2. Recording Time-Integral Current 
Density Distribution Using a 1-mn- 
Diameter-Aperture Faraday Cup 


This measurement was done on the 
short-pulse Febetron 706 system, 
whose cathode is circular. Figure 4 
shows its relative current density 
distribution. From this current 
density distribution, one can see 
some small nonuniformity. The cur- 
rent density at the center of the 
cathode is smaller than that at the 
edge. The peaks on both sides cor- 
respond to the 9-mm diameter which 
serves aS a measurement of the 
effective cathode size. 


| 


Relative current density 
~< 








9mm 














Radius 
Figure 4. Relative Current Density 
Profile Obtained by Smaller-—Aperture 
Faraday Cup 
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They are essentially 





3. Scanning Electron Micrographs of the Velvet Fibers After 1,000 Runs, Shown 
in Figure 5 [Photograph not reproduced] 


From the photographs, it is apparent that the tip of the fiber has melted. Its 
surface is relatively smooth and has been stretched. Another unique feature is 
that the velvet fibers become brittle after the experiment. They break off 
easily because they have been carbonized. However, they have not yet become 
conductors. 


III. Discussion of Results 

There are two theories that might be used to interpret the emission mechanism 
of the velvet cathode. One is that it is similar to an explosive emissior 
cathode. We believe that it might explain the mechanism in our experiment with 
longer pulses (>50 ns) because the electric field strength on the velvet is 
below 10°V/cem, which is the low-field case. Hence, electrons are initial] 
emitted from the tip of the fiber. Since the velvet is an insulator, the tip 
becomes charged after emission. 





For an explosive cathode, once the 25 
applied voltage exceeds the break- ’ @ Graphite ’ 


down voltage, the tips begin to 20F |e e Lead 
E I 
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the values are shown in Figure 6." Figure 6. Experimental Breakdown Dela 
The breakdown delay time is the Time as a Function of Appli¢« 
time after a pulse voltage is ap- Field and Electrode Material 
plied to a planar electrode to the 

formation of a flash point. Recently, we measured the breakdown delay tim f 
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Because of the huge difference in resistivity between metal and insulat 
the order of a dozen orders of magnitude, t, for velvet is very short This is 
the reason why the velvet cathode has avery short “closure time 


CL 





Another unique feature of the velvet cathode is that its plasma moves at a 
very slow speed. Theoretically, the velocity of the metal-—cathode plasma can 
be treated as a gas kinetics problem dealing with adiabatic diffusion of gas 
into vacuum. The velocity can be expressed as follows: 


v = (([4y/(y - 1)]) x n x E)*/? (2) 


where y = Cy. Cy which is the adiabatic coefficient, E is the latent heat of 
sublimation for the metal, and n is the overheat coefficient. Since the velvet 
cathode plasma temperature is lower than that of a metal cathode, the velvet 
cathode plasma moves at a much lower velocity than that of a metal cathode. 


How are electrons replenished from the metal substrate to the tips of the 
velvet cathode? There is a surface internal network theory. First, there are 
field-effect electrons in the high-field region near the metal—vacuum 
insulator interface. Then, these primary electrons bombard the surface of the 
insulator to create secondary electrons. Under the influence of the external 
field, electrons hop along the surface of the insulator to form a network on 
the velvet surface. Electrons are replenished to the tips by this mechanism. 
A detailed description of this mechanism is available in reference [9]. 


Recently, R.V. Latham, et al., presented an FIHEE (field-induced hot electron 
emission) model®’’:® to explain the prebreakdown electron emission mechanism on 
the high-voltage vacuum—insulator electrode surface. We believe that it might 
be used to explain the emission mechanism for the velvet cathode with short 
pulses (3 ns). 


With an applied field between the anode and cathode, velvet is an excellent 
dielectric material. Hence, the external field could penetrate the velvet and 
reach the metal substrate to cause the energy levels to change. 


There is a critical field in an insulator, approximately 10°~10°V/cm. When the 
applied field exceeds this strength, electrons can tunnel from the metal 
substrate to the conduction band of the velvet. This heats the velvet up and 
accelerates the ionization process which causes an avalanche in the conduction 
band of the velvet. 


Let us assume that AV is the voltage drop across the velvet fiber, V is the 
applied voltage between the anode and cathode, d is the distance between the 
two electrodes, Ad is the thickness of the velvet and ¢« is the relative 
dielectric constant of the velvet, then 


Ave (4) (42)v (3) 
€ 


In our experiment, ¢€ = 4, V = 300~600 kV, d = 8 mm, Ad = 1.5 mm. Here, AV = 
15~30 kV, and E = 10°~2 x 10°V/cm. Under such conditions, dielectric materials 
such as velvet have some conductive paths to convert the velvet from a 
nonconductive to a conductive state. 


5/7 





In the presence of a high field, the Poolo-Frenkel effect plays a major role 
(i.e., the field enhances thermal ionization between donors and traps). Then, 
the overall conductivity is a function of the applied field. In this case, J 
and AV may be expressed as follows: 


J 


i 


J,exp [ (e/kT) (ed V/ee,d) */?] (4) 


or J = J,exp[(e/2kT) (ed V/ee,d) '/*) 


depending upon the location and characteristics of the trap. 

















a 
In the short-pulse experiment, 
voltage and current profiles were 
measured after 30, 80, and 1/70 ae 
runs. From equation (3), we have 
the (AV)?/* versus lnI curve (Figure _ 
7). The results are in good agree- §& 
ment with the theory described 6+ ©No. |, after 80 shou 
above. No.2, after #0 shots 
GNo. 3, afler 170 shots 
In this case, the current density P P LL. } 
across the velvet is very high +130 130 140 130 160 170 


(approximately 7 kA/cm*); there- 
fore, some fibers are burnt. How- 
ever, the entire pulse lasts a very 
short time and burning might have 
occurred at the late stage of emis- 
sion. 


(AV) 
Figure 7. The Experimental Character- 
istics of lnI Versus (AV)?!/4 
of Velvet Cathode 
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(Text] Abstract 


The physical mechanism of microwaves excited by two relativistic electron 
beams (REBs) in a cylindrical waveguide is studied based on magnetic-—fluid 
theory. A dispersion relation has been derived in the presence of an infinite 
axial guiding magnetic field. An expression for the growth has also been found 
at the Cerenkov resonant point,’ i.e., w = kv,,.. The dependence of growth rate 
upon related parameters is also briefly discussed 


I. Introduction 


There are several types of high-power microwave generators,“ such as the 
dielectric Cerenkov microwave generator, virtual cathode oscillator (vircator) 
microwave generator, and beam—plasma interaction microwave generator. Each 
type of microwave generator has its own unique advantages and limitations in 
terms of band, power, and structure. 


A novel microwave generator, i.e., a dual-REB microwave generator, is 
presented in this paper. Physically speaking, when two electron beams move 
toward each other, they produce a two-stream instability.’ If the two beams are 


placed in a cylindrical waveguide, this instability can be coupled to yield 
waveguide electromagnetic modes to consequently generate a microwave output. 
Since the parameters of the two electron beams can be adjusted, it is easy to 


achieve resonance. 








II. Basic Equations 


Assume the initial velocities of the two cold REBs are vyy = ve, and Vg = 
—Vgo@z, respectively. The two beams are distinguished by subscripts a and #. 
The continuity equation is: 


ba + V+ (nyv,) = 0, (i=a, B) (1) 


The conservation of momentum equation is 


aP. 1 . 
J =- o <—- , = ) 2 
re e(—+ —v, xB), (i-=a, B) (2) 





where P; = myjvj represents the momentum of the electron; m, vj, and -e are the 
electron mass, velocity, and charge, respectively; c is the speed of light in 
vacuum; nj is the electron density, and yj = (1 - vi/c?)-}/? is the relativistic 
factor of the electron. E and B represent the electric and magnetic fields, 
respectively, and satisfy the Maxwell wave equation 


(w- 2 Sy gpe Vey ew (3) 


where J = » (-evjn;) is the total current density. 


Equations (1), (2), and (3) are the basic equations of the magnetic- fluid 
model. 


III. Dispersion Relation 


When an infinitely large axial guiding magnetic field is applied to a 
cylindrical waveguide, the two REBs only move in one-dimension along e,. Let 
us assume that the two REBs have the same radius r,, and the radius of the 
waveguide is R, the mean densities of the two beams are n,. and mgg, respec- 
tively. The dispersion relation may be derived from perturbation. To this end, 
we assume that 


Vi = Vio + Viz, Ni = Njo + nj,x, E = E,, B = B, (4) 
and . 
(Vii, Nii» E,, B,) ax e' (kz-vt) (5) 


where the subscript "1" represents a first-order perturbation term, and w and 
k are the angular frequency of the electromagnetic wave and wave number along 
e,, respectively. 





The two REBs in one-dimensional motion inside a cylindrical waveguide can only 
interact with the TM), mode, under the following boundary conditions: 


The following dispersion equation can be derived from equations (1)~(3): 


Y, (DR) Jg (pr,) - Jo(PR)¥g(pr,) — kK, dé (k,r,) 














._— (7) 
Y, (PR) J, (pr,) - J, (PR) Y, (pr,) Pp J,(kK,r;) 
where 
w? w? 
ki = [1 - 3 = - 7 p ] p? (8) 
Yao(KVgo-@)* = Yo (KVg.+@) ? 
2 
ps (9) 
c 


Jp and Y)y represent the zeroth-order Bessel function of the first and second 
kind, respectively. In equation (8) 


4ne*n., Vio 
pg © (ante, aig? (2 - 


IV. Growth Rate of Instability 


The dispersion equation (7) may be solved numerically. Let A,,/R be the nth 
root (n=1,2---) of k, in equation (7). When R = rp, Ap, is the root of Jog(Agn) 
= 0. Hence, from equation (8), we have 











2 2 os 
(1 - oe ° Ory (uy? - K2c2) = Ame (10) 
Yeo (@-kv,,)? Ypo (@+KvV,,)? R? 


Equation (10) clearly points out that the physical process in a dual—REB 
microwave generator is the process that leads to the coupling of the 
electromagnetic mode of the waveguide due to two-stream instability. 


If both beams have similar electron velocities and densities, it is possible 
to find the growth rate of the pure growth mode where the real part of the 
frequency is close to zero from equation (10). Since the frequency is near 
zero, this mode is of no interest to us. Hence, let us consider the following 
situation where the density of electron beam a is much smaller than that of 
electron beam §, i.e., 








2 2 
a « = (11) 
Yao YBo 


Near the Cerenkov resonant point w = Kvy), because there might be some strong 
interaction between the electromagnetic wave and electron beama, the maximum 
growth rate of instability can be obtained from equation (10). 
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+ 1/3 
Wie | Veo Vp0 ) 
- v3 Veo ; 
— 24Bye | 12 6c? V, 7 te) 
e8 1 wo+ —— (yao-1) [ (yeo-2) (1+ P2) +1] 
L R Veo 
where 
: 1/2 
_ Ws _ donYao Yao 
™ ° 3 Vao R*4 (13) 
Ypo(1+—F°)? 
«0 
The frequency of the microwave is 
wo = wo, - T/¥3 = w, (14) 


V. Conclusions 


A magnetic-fluid model is used to study the interaction between two electron 
beams and a waveguide. Under the approximation that the densities of the two 
electron beams differ significantly, an expression, equation (12), for the 
maximum growth rate near the Cerenkov resonant point, w = kv,,, is obtained. 
The results indicate that: 


(1) The output microwave frequency is a function of vVy4o, Vgo, Mao, Ngo, R and 

rh. By properly adjusting some parameters while keeping others unchanged, it 
is possible to achieve some degree of resonance. It might be more effective to 
adjust the electron density ng, of REB 6 and the electron velocity v,, of REB 


a. 
(2) The growth rate can be increased by raising the n,,. and ng of the REBs. 


(3) The growth rate can be lowered by increasing the velocity, or energy, of 
the REBs. 


2 2 2 2 
W op AonYe0Vao 


tii «7. ) <is een ; 
, Veo R2 when n increases, Ap, goes up and the growth 
(1+—~)? 
- 





(4) With 
Yo 

@\ 
rate of instability falls. Therefore, the growth rate of the lowest-order mode 


TM,, should be the highest. 


(5) In order to illustrate that the electromagnetic radiation emitted by the 
dual-REB microwave generator is in the microwave range, a numerical example is 
given below. 


Assume REB a and REB § have the same velocity, i.e., Yao = Yeo = 1.196 and vy 
= Veo = 0.548c. However, the densities are quite different: n,. = 2.688 x 
10'4em™? and ng = 4.30 x 10'4cm™’. The radius of the cylindrical waveguide 

R= 1.5 cm. The electron beams almost fill the waveguide (r,~R). From eq. (1 
we can calculate that w, = 3.178 x 10’°s"'. From w. = (2”/A)v,,, the wavelength 
of the electromagnetic wave was calculated to be 3.248 cm. 
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